d— wave collapse and explosion of a dipolar Bose-Einstein condensate 
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We investigate the collapse dynamics of a dipolar condensate of 52 Cr atoms when the s— wave 
scattering length characterizing the contact interaction is reduced below a critical value. A complex 
dynamics, involving an anisotropic, d— wave symmetric explosion of the condensate, is observed. 
The atom number decreases abruptly during the collapse. We find good agreement between our 
experimental results and a those of a numerical simulation of the three-dimensional Gross-Pitaevskii 
equation, including contact and dipolar interactions as well as three-body losses. The simulation 
indicates that the collapse induces the formation of two vortex rings with opposite circulations. 
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The underlying symmetries of physical systems often 
determine the nature and dynamics of macroscopic quan- 
tum states. For example, the difference between isotropic 
and d— wave pairing of electrons, in conventional and 
high-T c superconductors, respectively, leads to funda- 
mentally different properties [1]. Degenerate quantum 
gases are usually dominated by isotropic (s— wave) con- 
tact interactions. Dipolar quantum gases (i.e. in which 
the dipole-dipole interaction (DDI) between permanent 
dipole moments play a significant or even dominant role) 
are governed by the d— wave symmetry of the long-range 
DDI, which gives rise to novel properties. 

Examples of fascinating predictions for polarized dipo- 
lar quantum gases range from a roton-maxon spec- 
trum [2] for the elementary excitations of a quasi two- 
dimensional dipolar Bose-Einstein condensate (BEC), to 
the existence of novel quantum phases (such as a 'checker- 
board' insulator, or a supersolid) for dipolar quantum 
gases in optical lattices [3| . The DDI also modifies the hy- 
drodynamic equations describing the dynamics of a BEC, 
which has been probed experimentally by studying the 
expansion of the cloud when released from the trap [4|. 
UnusuaL structured shapes for the BEC have been pre- 
dicted [1,0. In the unpolarized case, the DDI dramati- 
cally enriches the physics of spinor BECs [7]. Ultracold 
dipolar fermions also have fascinating properties 

A striking example of the new properties of dipolar 
BECs is given by their stability, which, contrary to the 
case of contact interaction, depends strongly on the trap 
geometry. Consider a pancake-shaped trap with the 
dipole moments of the particles oriented perpendicular 
to the plane of the trap. The DDI is then essentially re- 
pulsive, and the BEC is stable, independently of the atom 
number. In contrast, a cigar-shaped trap cannot stabi- 
lize a purely dipolar BEC. We experimentally studied @ 
this geometry-dependent stability of a dipolar quantum 
gas by using a 52 Cr BEC, and mapped out the stability 



diagram of the condensate as a function of the scatter- 
ing length a (characterizing the contact interaction) and 
the trap aspect ratio. In the case of a pure contact in- 
teraction, crossing the stability border into the unstable 
regime a < leads to a collapse of the BEC [Io|, EH, EH. 
This gives rise to an interesting dynamics involving a fast 
implosion of the condensate followed by the formation of 
energetic 'bursts' of atoms [13], or the formation of soli- 
ton trains 0,0. 

In this Letter, we investigate experimentally the col- 
lapse dynamics of a dipolar 52 Cr BEC when the scat- 
tering length a is decreased (by means of a Feshbach 
resonance) below the critical value a cr i t for stability 0. 
We observe a rich dynamics on a timescale shorter than 
the trap period, with the formation of an expanding 
structure featuring a d— wave symmetry. We study the 
atom number in the condensate as a function of time, 
and find an abrupt decrease due to inelastic losses. Fi- 
nally, we compare our experimental results with a three- 
dimensional numerical simulation of the Gross-Pitaevskii 
equation (GPE) including both contact and dipolar in- 
teractions, as well as three-body losses. Such a general- 
ized GPE with three-body losses has been demonstrated 
to explain the main features of 'Bose-nova' experiments 
with 85 Rb [lfl]. Here we generalize this model to include 
the DDI. As shown later, this generalized model accounts 
very well for the observed d— wave collapse. 

The experimental setup to produce a 52 Cr BEC above 
the Feshbach resonance located at a magnetic field Bq ~ 
589 G has been described elsewhere [9]. Close to the 
resonance, the scattering length a varies with the applied 
magnetic field B as 



a(B) = a bg 1 



B-B 



where A ~ 1.5 G is the resonance width, and a bg — 
100 ao the background scattering length (ao is the Bohr 
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FIG. 1: (color). Collapse dynamics of the dipolar condensate, 
(a) Timing of the experiment. The red curve represents the 
time variation of the scattering length a(t) one would have in 
the absence of eddy currents, while the blue curve is obtained 
by taking them into account (see text), (b) Sample absorption 
image of the collapsed condensate for thoid = 0.4 ms, after 
8 ms of time of flight, showing a 'clover leaf pattern on top of 
a broad thermal cloud. This image was obtained by averaging 
60 pictures taken under the same conditions, (c) Same image 
as (b) with the thermal cloud subtracted. In (b) and (c) the 
field of view is 270 /im by 270 /im. The green arrow indicates 
the direction of the magnetic field, (d) Series of images of the 
condensate for different values of thoid (upper row) and results 
of the numerical simulation without adjustable parameters 
(lower row); the field of view is 130 /im by 130 /im. 



radius). We calibrate the variation a(B) of the scattering 
length by measuring, after expansion, the BEC size and 
atom number [9[ . The reduction of a close to Bo + A is 
accompanied by inelastic losses. By measuring the 1/e 
lifetime and the density of the BEC close to resonance, we 
estimate the three-body loss coefficient to be constant for 
the range of scattering lengths (5 ^ a/ao ^ 30) studied 
here, with a value L 3 ~ 2 x 10 -40 m 6 /s. 

To study the collapse dynamics, we first create a 
BEC of typically 20,000 atoms in a trap with frequen- 
cies (v x ,v y ,v z ) ~ (660,400,530) Hz at a magnetic field 
~ 10 G above the Feshbach resonance, where the scatter- 
ing length is a ~ 0.9 a Dg . We then decrease a by ramping 
down B linearly over 8 ms to a value cl[ = 30 ao which still 
lies well above the critical value for collapse, measured to 
be at a cr it — (15 ±3) ao [shaded area on Fig. 1(a)] for our 
parameters [9[ . This ramp is slow enough to be adiabatic 
(a I a <C v x ,y,z), so that the BEC is not excited during it. 
After 1 ms waiting time, a is finally ramped down rapidly 
to <2f = 5ao, which is below the collapse threshold. For 
this, we ramp linearly in 1 ms the current I(t) in the coils 
providing the magnetic field B = al. However, due to 
eddy currents in the metallic vacuum chamber, the ac- 
tual value of B(t) and hence that of a(t) change in time 
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FIG. 2: (color). In-trap column density obtained in the sim- 
ulation, for different thoid- The field of view is 5 /im by 5 
/im. Due to the DDI, the condensate collapses radially, ac- 
quiring the shape of a very thin cigar elongated along z. At 
thoid — 0.5 ms, the collapse occurs, and immediately after, 
the cloud starts to expand radially. 

as depicted in blue on Fig. 1(a). To obtain this curve, we 
used Zeeman spectroscopy to measure the step response 
of B(t) to a jump in the current I(t) (corresponding to a 
~ 15 G change in B), and found that the resulting B(t) is 
well described if tB + B = al(t) holds, with r ~ 0.5 ms. 
From this equation and the measured I(t) we determine 
the actual a(t). 

After the ramp, we let the system evolve for an ad- 
justable time thoid an d then the trap is switched off. Note 
that the origin of thoid corresponds to the end of the 
ramp in I(t). Because of eddy currents, thoid = about 
0.2 ms before the time at which the scattering length 
crosses a cr it- However, as we shall see below, even for 
^hoid < 0.2 ms a collapse (happening not in trap, but dur- 
ing the time of flight) is observed, since during expansion 
the scattering length continues to evolve towards <2f . The 
large magnetic field along z is rapidly turned off (in less 
than 300 /is) after 4 ms of expansion, and the conden- 
sate expands for another 4 ms in an 11 G field pointing 
in the x direction, before being imaged by absorption of 
a resonant laser beam propagating along x. Changing 
the direction of the field allows us to use the maximum 
absorption cross-section for the imaging (if the latter was 
done in high field, the absorption cross-section would be 
smaller, thus reducing the signal to noise ratio of the 
images). We checked that this fast switching has no in- 
fluence on the condensate dynamics. We observe that 
the atomic cloud has a clear bimodal structure, with a 
broad isotropic thermal cloud, well fitted by a Gaussian, 
and a much narrower, highly anisotropic central feature, 
interpreted as the remnant BEC [see Fig. 1(b) and (c)]. 

The upper row of Fig. 1(d) shows the time evolu- 
tion of the condensate when varying thoid- The images 
were obtained by averaging typically five absorption im- 
ages taken under the same conditions; the thermal back- 
ground was subtracted, and the color scale was adjusted 
separately for each thoid for a better contrast. From an 
initial shape elongated along the magnetization direc- 
tion z, the condensate rapidly develops a complicated 
structure with an expanding, torus-shaped part close to 
the z = plane. Interestingly, the angular symme- 
try of the condensate at some specific times (e.g. at 
^hoid = 0.5 ms) is reminiscent of the d— wave angular 
symmetry 1 — 3 cos 2 9 of the DDI. For larger values of 
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^hoid? we observe that the condensate 'refocuses' due to 
the presence of the trap [l7| . 

The lower row of Fig. 1(d) shows the column density 
f \ip(r)\ 2 dx (where ip(r) is the order parameter of the 
condensate after time of flight) obtained from a numerical 
simulation of the three-dimensional GPE 
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where 



[/(r,t) = ^ (r) 



fio/a 2 1 — 3 cos 2 
4tt r 3 



stands for the contact and dipolar interactions, 6 being 
the angle between r and the direction of polarization. 
Here m is the atomic mass, fio the permeability of vac- 
uum, and fi = 6/iB the magnetic moment of a Cr atom 
(/iB is the Bohr magneton). The non-unitary term pro- 
portional to Ls describes three-body losses. The scatter- 
ing length a(t) is changed according to the blue curve in 
Fig. 1(a) and the trap potential Vt rap is switched off 
at the beginning of the 8 ms time of flight. For the 
simulation, space is discretized into a 128 x 128 x 128 
mesh with a step size of 70 nm. For the kinetic part, 
the Crank-Nicolson scheme is used for the time evolution 
to avoid numerical instability. For the interaction part, 
the convolution integral is calculated using a fast Fourier 
transform. After the trap is switched off, the mesh is 
extended to 512 x 512 x 512 to describe the expansion 
of the cloud. When the density becomes low enough so 
that the nonlinear terms of the GPE can be neglected, 
the free expansion propagator is used to give the final 
time-of-flight images. The agreement between the exper- 
imental data and the simulation, performed without any 
adjustable parameter, is excellent. 

The cloverleaf patterns seen in Fig. 1(b)- (d) are caused 
by the anisotropic collapse and the subsequent dynamics 
of the system. Figure 2 shows the in-trap evolution of 
the condensate as a function of thoid- The mechanism 
of the condensate 'explosion' is as follows [13]: When 
the atomic density grows due to the attractive interac- 
tion, three-body losses predominantly occur in the high- 
density region. The centripetal force is then decreased, 
and the atoms that gathered in this narrow central re- 
gion are ejected due to the 'quantum pressure' arising 
from the uncertainty principle. The kinetic energy is 
supplied by the loss of the negative interaction energy. 
For the contact interaction, the collapse and subsequent 
atomic 'explosion' is isotropic [18]. In the present case, 
the collapse occurs mainly in the x — y direction due to 
anisotropy of the DDI (in the absence of inelastic losses, 
the condensate would indeed become an infinitely thin 
cigar-shaped cloud along z, see Fig. 3 of Ref.[9], and the 
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FIG. 3: (color). Atom losses during collapse. Blue cir- 
cles: atom number A^bec in the condensate as a function 
of thoid- The solid curve is the result of the simulation for 
L3 = 2 x 10 -40 m 6 /s, without any adjustable parameter. 
Inset: the atom number iVthermai in the thermal cloud (red 
squares) remains essentially constant during the collapse. 



in-trap image at thoid = 0.5 ms in Fig. 2), and therefore 
the condensate 'explodes' essentially radially, producing 
the anisotropic shape of the cloud. The numerical sim- 
ulation reveals that, for thoid < 0.5 ms, the collapse ob- 
served in Fig. 1(d) occurs not during the holding time 
but during the time-of-flight. We stress that in the ab- 
sence of three-body losses, the explosion following the 
collapse would not be observed. 

From the images, the atom number A/bec in the con- 
densate is obtained by integrating the optical density. 
Blue circles in Fig. 3 show TVbec as a function of thoid- 
The BEC atom number is initially 7V B ec(0) ~ 16,000 
and decreases toward its asymptotic value ~ 6, 000. Over 
the same time scale, the atom number iVthermai in the 
thermal cloud (inset of Fig. 3) stays constant. The size 
of the thermal cloud after expansion is also constant over 
this period. This suggests that the thermal cloud does 
not play any significant role in the collapse dynamics. 
For thoid < 0.5 ms, the collapse actually occurs during 
the time of flight, which explains the gradual decay of 
^BEc(^hoid), and why the atom losses are not as large 
as those when the collapse occurs in trap (in the latter 
case, 70% of the atoms are lost). The missing atoms have 
very likely escaped from the trap as energetic molecules 
and atoms produced in three-body collisions. This is con- 
firmed by the fact that the simulation gives a A^BEc(^hoid) 
curve (solid line in Fig. 3) which matches well the exper- 
imental data. Experimental uncertainties in the parame- 
ters used in the simulation (trap frequencies, values of L3 
and r) probably explain the small discrepancy between 
the experiment and the numerical results. 

The numerical simulation gives access not only to the 
density |^(r)| 2 , but also to the phase S(r) of the order 
parameter ip (i.e. to the velocity field v = KVS/m) and 
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FIG. 4: (color). Vortex rings predicted by the numerical sim- 
ulation, (a) Iso-density surface of an in-trap condensate at 
thoid = 0.8 ms. The topological defects are shown by the red 
rings, (b) Velocity field of the atomic flow in the x = plane 
at thoid — 0-8 m s. The field of view is 2.5 /mi by 2.5 /mi; the 
color scale represents the velocity (red is faster). 

reveals the generation of vortex rings [H, Q of charge 
±1. Figure 4(a) shows an in-trap iso-density surface of 
a condensate at thoid = 0.8 ms and the location of the 
vortex rings (shown as red curves). Comparing Fig. 4(a) 
and Fig. 2 (at thoid = 0.8 ms )? we find that the topologi- 
cal defects encircle the two 'leaves' of the 'clover'. Figure 
4(b) shows the velocity field v(r) in the x = plane. The 
atoms ejected in the x — y plane flow outward, while the 
atoms near the z axis still flow inward, giving rise to the 
circulation. Thus, the vortex-ring formation is specific to 
the d— wave collapse induced by the DDL Although the 
vortex rings are not observed directly in the experiment, 
the excellent agreement between the experiment and sim- 
ulation in Figs. 1 and 3 strongly suggests the creation of 
vortex rings during the collapse of the condensate. 

In conclusion, we have investigated the collapse dy- 
namics of a dipolar BEC. Contrary to the case of an 
isotropic contact interaction, the DDI induces the for- 
mation of a structured cloud featuring a d— wave sym- 
metry. The collapse dynamics is quantitatively repro- 
duced by numerical simulations of the GPE without any 
adjustable parameter. An interesting subject for future 
studies is the dependence of the collapse dynamics on 
the trap geometry: one may wonder if the condensate 
would collapse in the same way if initially trapped in a 
very anisotropic (e.g., pancake-shaped) trap. A natural 
extension of this work would involve detecting, e.g. by 
interferometric methods [U [22j], the vortex rings pre- 
dicted by the simulation. Finally, whether one can nu- 
cleate stable vortex rings by initiating the collapse, and 
then changing a back to a value corresponding to a sta- 
ble BEC, is a question which certainly deserves further 
investigations. 
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